Search for: All records

ABSTRACT The dark matter (DM) distribution in dwarf galaxies provides crucial insights into both structure formation and the particle nature of DM. GraphNPE (Graph Neural Posterior Estimator), first introduced in Nguyen et al. (2023), is a novel simulation-based inference framework that combines graph neural networks and normalizing flows to infer the DM density profile from line-of-sight stellar velocities. Here, we apply GraphNPE to satellite dwarf galaxies in the FIRE-2 Latte simulation suite of Milky Way-mass haloes, testing it against both Cold and Self-Interacting DM scenarios. Our method demonstrates superior precision compared to conventional Jeans-based approaches, recovering DM density profiles to within the 95 per cent confidence level even in systems with as few as 30 tracers. Moreover, we present the first evaluation of mass modelling methods in constraining two key parameters from realistic simulations: the peak circular velocity, $$V_\mathrm{max}$$, and the peak virial mass, $$M_\mathrm{200m}^\mathrm{peak}$$. Using only line-of-sight velocities, GraphNPE can reliably recover both $$V_\mathrm{max}$$ and $$M_\mathrm{200m}^\mathrm{peak}$$ within our quoted uncertainties, including those experiencing tidal effects ($$\gtrsim 63~{{\rm per\ cent}}$$ of systems are recovered within our 68 per cent confidence intervals and $$\gtrsim 92~{{\rm per\ cent}}$$ within our 95 per cent confidence intervals). The method achieves $$10-20~{{\rm per\ cent}}$$ accuracy in $$V_\mathrm{max}$$ recovery, while $$M_\mathrm{200m}^\mathrm{peak}$$ is recovered to $$0.1-0.4 \, \mathrm{dex}$$ accuracy. This work establishes GraphNPE as a robust tool for inferring DM density profiles in dwarf galaxies, offering promising avenues for constraining DM models. The framework’s potential extends beyond this study, as it can be adapted to non-spherical and disequilibrium models, showcasing the broader utility of simulation-based inference and graph-based learning in astrophysics. 

Abstract Spatial patterns of stellar elemental abundances encode rich information about a galaxy’s formation history. We analyze the radial, vertical, and azimuthal variations of metals in stars, both today and at formation, in the FIRE-2 cosmological simulations of Milky Way (MW)-mass galaxies, and we compare them with the MW. The radial gradient today is steeper (more negative) for younger stars, which agrees with the MW, although radial gradients are shallower in FIRE-2. Importantly, this age dependence was present already at birth: radial gradients today are only modestly (≲0.01 dex kpc⁻¹) shallower than at birth. Disk vertical settling gives rise to negative vertical gradients across all stars, but vertical gradients of mono-age stellar populations are weak. Similar to the MW, vertical gradients in FIRE-2 are shallower at larger radii, but they are overall shallower in FIRE-2. This vertical dependence was present already at birth: vertical gradients today are only modestly (≲0.1 dex kpc⁻¹) shallower than at birth. Azimuthal scatter is nearly constant with radius, and it is nearly constant with age ≲8 Gyr ago but increases for older stars. Azimuthal scatter is slightly larger (≲0.04 dex) today than at formation. Galaxies with larger azimuthal scatter have a stronger radial gradient, implying that azimuthal scatter today arises primarily from the radial redistribution of gas and stars. Overall, spatial variations of stellar metallicities show only modest differences between formation and today; spatial variations today primarily reflect the conditions of stars at birth, with spatial redistribution of stars after birth contributing secondarily. 

Abstract Using the FIRE-2 cosmological zoom-in simulations, we investigate the temporal evolution of gas-phase metallicity radial gradients of Milky Way–mass progenitors in the redshift range of 0.4 <z< 3. We pay special attention to the occurrence of positive (i.e., inverted) metallicity gradients—where metallicity increases with galactocentric radius. This trend, contrary to the more commonly observed negative radial gradients, has been frequently seen in recent spatially resolved grism observations. The rate of occurrence of positive gradients in FIRE-2 is about ∼7% for 0.4 <z< 3 and ∼13% at higher redshifts (1.5 <z< 3), broadly consistent with observations. Moreover, we investigate the correlations among galaxy metallicity gradient, stellar mass, star formation rate (SFR), and degree of rotational support. Metallicity gradients show a strong correlation with both sSFR and the rotational-to-dispersion velocity ratio (v_c/σ), implying that starbursts and kinematic morphology of galaxies play significant roles in shaping these gradients. The FIRE-2 simulations indicate that galaxies with high sSFR ( $\log \left(\mathrm{sSFR}\right[{\mathrm{yr}}^{-1}\left]\right)\gtrsim -9.2$) and weak rotational support (v_c/σ≲ 1) are more likely—by ∼15%—to develop positive metallicity gradients. This trend is attributed to galaxy-scale gas flows driven by stellar feedback, which effectively redistribute metals within the interstellar medium. Our results support the important role of stellar feedback in governing the chemo-structural evolution and disk formation of Milky Way–mass galaxies at the cosmic noon epoch. 

ABSTRACT We analyse the stellar distributions on the [Fe/H]–[Mg/Fe] plane for 11 Milky Way-mass galaxies from the FIRE-2 cosmological baryonic zoom-in simulations. Alpha-element bimodality, in the form of two separate sequences on the [Fe/H]–[Mg/Fe] plane, is not a universal feature of disc galaxies. Five galaxies demonstrate double sequences with the $$\alpha$$-enriched one being older and kinematically hotter, in qualitative agreement with the high-$$\alpha$$ and low-$$\alpha$$ populations in the Milky Way disc; three galaxies have unimodal distribution, two show weakly bimodal features where low-$$\alpha$$ sequence is visible only over a short range of metallicities, and one show strong bimodality with a different slope of high-$$\alpha$$ population. We examine the galaxies’ gas accretion history over the last 8 Gyr, when bimodal sequences emerge, and demonstrate that the presence of the low-$$\alpha$$ sequence in the bimodal galaxies is related to the recent infall of metal-poor gas from the circumgalactic medium that joins the galaxy in the outskirts and induces significant growth of the gas discs compared to their non-bimodal counterparts. We also analyse the sources of the accreted gas and illustrate that both gas-rich mergers and smooth accretion of ambient gas can be the source of the accreted gas, and create slightly different bimodal patterns. 

Abstract The physical mechanisms responsible for bar formation and destruction in galaxies remain a subject of debate. While we have gained valuable insight into how bars form and evolve from isolated idealized simulations, in the cosmological domain, galactic bars evolve in complex environments, with mergers and gas accretion events occurring in the presence of the turbulent interstellar medium with multiple star formation episodes, in addition to coupling with their host galaxies’ dark matter halos. We investigate the bar formation in 13 Milky Way–mass galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulations. 8 of the 13 simulated galaxies form bars at some point during their history: three from tidal interactions and five from internal evolution of the disk. The bars in FIRE-2 are generally shorter than the corotation radius (mean bar radius ∼1.53 kpc), have a wide range of pattern speeds (36–97 km s⁻¹kpc⁻¹), and live for a wide range of dynamical times (2–160 bar rotations). We find that the bar formation in FIRE-2 galaxies is influenced by satellite interactions and the stellar-to-dark-matter mass ratio in the inner galaxy, but neither is a sufficient condition for bar formation. Bar formation is more likely to occur, with the bars formed being stronger and longer-lived, if the disks are kinematically cold; galaxies with high central gas fractions and/or vigorous star formation, on the other hand, tend to form weaker bars. In the case of the FIRE-2 galaxies, these properties combine to produce ellipsoidal bars with strengthsA₂/A₀∼ 0.1–0.2. 

We construct time-evolving gravitational potential models for a Milky Way–mass galaxy from the FIRE-2 suite of cosmological-baryonic simulations using basis function expansions. These models capture the angular variation with spherical harmonics for the halo and azimuthal harmonics for the disk, and the radial or meridional plane variation with splines. We fit low-order expansions (four angular/harmonic terms) to the galaxy’s potential for each snapshot, spaced roughly 25 Myr apart, over the last 4 Gyr of its evolution, then extract the forces at discrete times and interpolate them between adjacent snapshots for forward orbit integration. Our method reconstructs the forces felt by simulation particles with high fidelity, with 95% of both stars and dark matter, outside of self-gravitating subhalos, exhibiting errors ≤4% in both the disk and the halo. Imposing symmetry on the model systematically increases these errors, particularly for disk particles, which show greater sensitivity to imposed symmetries. The majority of orbits recovered using the models exhibit positional errors ≤10% for 2–3 orbital periods, with higher errors for orbits that spend more time near the galactic center. Approximate integrals of motion are retrieved with high accuracy even with a larger potential sampling interval of 200 Myr. After 4 Gyr of integration, 43% and 70% of orbits have total energy and angular momentum errors within 10%, respectively. Consequently, there is higher reliability in orbital shape parameters such as pericenters and apocenters, with errors ∼10% even after multiple orbital periods. These techniques have diverse applications, including studying satellite disruption in cosmological contexts. 

Abstract Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, pulsar timing, and eclipsing binary stars in the solar neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density and infer the amount of dark matter (DM) by accounting for the mass in stars, gas, and dust. Leveraging FIRE-2 simulations of Milky Way–mass galaxies we compare vertical acceleration profiles between cold DM (CDM) and self-interacting DM (SIDM) with a constant cross section of 1 cm²g⁻¹across three halos with diverse assembly histories. Notably, significant asymmetries in vertical acceleration profiles near the midplane at fixed radii are observed in both CDM and SIDM, particularly in halos recently affected by mergers with satellites of Sagittarius/SMC-like masses or greater. These asymmetries offer a unique window into exploring the merger history of a galaxy. We show that SIDM halos manifest a more oblate shape and consistently exhibit higher local stellar and DM densities and steeper vertical acceleration gradients, up to 10%–30% steeper near the solar neighborhood. However, similar magnitude changes can arise from azimuthal variations in the baryonic components at a fixed radius and external influences like mergers, making it difficult to distinguish between CDM and SIDM using acceleration measurements in a single galaxy. 

Abstract Open-star clusters are the essential building blocks of the Galactic disk; “strong chemical tagging”—the premise that all star clusters can be reconstructed given chemistry information alone—is a driving force behind many current and upcoming large Galactic spectroscopic surveys. In this work, we characterize the abundance patterns for nine elements (C, N, O, Ne, Mg, Si, S, Ca, and Fe) in open clusters (OCs) in three galaxies (m12i, m12f, and m12m) from the Latte suite of FIRE-2 simulations, to investigate the feasibility of strong chemical tagging in these simulations. We select young massive (≥10^4.6M_⊙) OCs formed in the last ∼100 Myr and calculate the intra- and intercluster abundance scatter for these clusters. We compare these results with analogous calculations drawn from observations of OCs in the Milky Way. We find the intracluster scatter of the observations and simulations to be comparable. While the abundance scatter within each cluster is minimal (≲0.020 dex), the mean abundance patterns of different clusters are not unique. We also calculate the chemical difference in intra- and intercluster star pairs and find it, in general, to be so small that it is difficult to distinguish between stars drawn from the same OC or from different OCs. Despite tracing three distinct nucleosynthetic families (core-collapse supernovae, white dwarf supernovae, and stellar winds), we conclude that these elemental abundances do not provide enough discriminating information to use strong chemical tagging for reliable OC membership. 

Abstract Stellar feedback influences the star formation rate (SFR) and the interstellar medium of galaxies in ways that are difficult to quantify numerically, because feedback is an essential ingredient of realistic simulations. To overcome this, we conduct a feedback-halting experiment starting with a Milky Way–mass galaxy in the second-generation Feedback In Realistic Environments (FIRE-2) simulation framework. By terminating feedback, and comparing to a simulation in which feedback is maintained, we monitor how the runs diverge. We find that without feedback, the interstellar turbulent velocities decay. There is a marked increase of dense material, while the SFR increases by over an order of magnitude. Importantly, this SFR boost is a factor of ∼15–20 larger than is accounted for by the increased freefall rate caused by higher densities. This implies that feedback moderates the star formation efficiency per freefall time more directly than simply through the density distribution. To probe changes at the scale of giant molecular clouds (GMCs), we identify GMCs using density and virial parameter thresholds, tracking clouds as the galaxy evolves. Halting feedback stimulates rapid changes, including a proliferation of new bound clouds, a decrease of turbulent support in loosely bound clouds, an overall increase in cloud densities, and a surge of internal star formation. Computing the cloud-integrated SFR using several theories of turbulence regulation, we show that these theories underpredict the surge in SFR by at least a factor of 3. We conclude that galactic star formation is essentially feedback regulated on scales that include GMCs, and that stellar feedback affects GMCs in multiple ways.